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ABSTRACT 

We use simulated maps of the cosmic microwave background anisotropy to quan- 
tify the ability of different statistical tests to discriminate between Gaussian and non- 
Gaussian models. Despite the central limit theorem on large angular scales, both the 
genus and extrema correlation are able to discriminate between Gaussian models and 
a semi-analytic texture model selected as a physically motivated non-Gaussian model. 
When run on the COBE 4-year CMB maps, both tests prefer the Gaussian model. 
Although the bispectrum has comparable statistical power when computed on the full 
sky, once a Galactic cut is imposed on the data the bispectrum loses the ability to 
discriminate between models. Off-diagonal elements of the bispectrum are comparable 
to the diagonal elements for the non-Gaussian texture model and must be included to 
obtain maximum statistical power. 

Subject headings: cosmic microwave background — cosmology: observations — meth- 
ods: statistical 



1. Introduction 

The statistics of the cosmic microwave background (CMB) probe physical conditions in the 
early universe. Most inflationary models predict the CMB to follow Gaussian statistics (Bardeen et. al. 1983 
Guth and Pi 1982; Hawking 1982; Starobinskii 1982). Searching CMB data for non-Gaussian dis- 
tributions is thus an important test of inflationary cosmologies. Recent analyses of the full-sky 
COBE data provide mixed results. Kogut et. al. (Kogut et. al. 1995; Kogut et. al. 1996) use 
a likelihood analysis to compare the COBE-DMR 2- and 4-year maps to the standard Gaussian 
model and a set of non-Gaussian toy models. The COBE data lie near the mode of the sta- 
tistical distributions of the genus, 3-point correlation, and correlation of extrema points for the 
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Gaussian model, while rejecting several non-Gaussian alternatives. Other authors using similar 
statistics on the COBE maps find comparable results (Hinshaw et. al. 1994; Smoot et. al. 1994; 
Luo 1994; Torres et al. 1995; Colley et. al. 1996; Heavens 1998). More recently, several authors 
have challenged the conclusion that the COBE data are well described by Gaussian statistics. 
Tests using the normalized bispectrum (Ferreira et. al. 1998; Magueijo 2000) and a wavelet anal- 
ysis (Pando et. al. 1998) both show the COBE data to be incompatible with Gaussian simula- 
tions at 98% confidence; however, these analyses did not test any non-Gaussian alternatives. 
The recent work of Mukherjee (2000) brings into question the the statistical significance of the 
(Pando et. al. 1998) result. In (Mukherjee et. al. 2000) a proper Monte-Carlo estimation of the 
significance level shows the wavelet technique does not yield strong evidence for non-Gaussianity in 
the 4-year COBE data. Bromley and Tegmark (1999) argue that the non-Gaussian signal (if any) 
must lie in the phase information in the COBE sky maps, while Banday et al. (2000) ascribe the 
discrepancy to instrumental artifacts in one of the DMR frequency channels. 

Searches for non-Gaussianity are complicated by several factors. On large angular scales, the 
superposition of multiple sources within the instrument angular resolution tends to push initially 
non-Gaussian distributions toward the Gaussian model (the Central Limit Theorem). Sample 
variance can occasionally produce a non-Gaussian distribution drawn from a Gaussian parent pop- 
ulation. Given enough independent tests of the same observable sky, some outliers are expected 
(Bromley and Tegmark 1999). Demonstrating that a statistic calculated for a CMB data set lying 
far from the mode expected for a Gaussian sky is a necessary but not sufficient condition to conclude 
that the data are not, in fact, a sample drawn from a Gaussian parent population. 

The statistical power of a test quantifies the probability of correctly accepting or rejecting 
hypotheses based on the results of that test. Quantifying the power for a given statistic thus 
requires at least two competing hypotheses. Since there is a single Gaussian distribution but an 
infinite number of non-Gaussian distributions, we must necessarily restrict the scope of specific 
non-Gaussian hypotheses. (Kogut et. al. 1995; Kogut et. al. 1996) used non-Gaussian toy models 
to show that the genus and extrema correlation could discriminate between Gaussian and non- 
Gaussian populations. We extend this analysis to include the bispectrum, replacing the earlier 
toy models with a model of topological defects as a specific example of a physically-motivated 
alternative to the standard Gaussian model. 

2. Likelihood Analysis and Statistical Power 

Maximum-likelihood techniques are a commonly used tool for statistical inference. Given any 
set of K statistics derivable from a sky map, which we view as the vector S, we define 

K 

X 2 (S) = ^ (S a - (S^iM-XpiSp - {Sp)) (1) 

a,/3=l 
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where (S) is the mean value of the statistic for the model under consideration, and M its covariance 
matrix. The likelihood is then 

C (S) = -^EtMl (2) 
^/(2tt) k det(M) 

where we make the convenient assumption that the statistics S obey a multivariate Gaussian 
probability distribution. Since the large-scale CMB anisotropy is close to a Gaussian distribution, 
We do not expect the calculation of the true distributions based on Monte-Carlo simulations can 
significantly change the relative statistical powers we determine. For the simple comparison of two 
models, we compute the likelihood functions L\ and £2 for a given statistic S(AT). The data AT 
are more likely to be drawn from parent population model 1 if £i(S(AT)) > £2(8 (AT)), and vice 
versa for model 2. 

For a statistic to be useful, we need to determine how well it can discriminate between the 
models under consideration. Recall there are two types of error: a Type I error rejects a hypothesis 
when it is correct, while a Type II error accepts a hypothesis when it is false. The statistical power 
of a test quantifies the robustness of a test against Type II errors. By sampling the likelihood 
functions £j with independent samples drawn from each of the two models under consideration, 
we obtain the probability to incorrectly identify a sample from model i as drawn from model j: 
P(Si\Hj), where Sj is extracted from a realization of model i. P(Si\Hi) is thus the probability 
to correctly identify a sample from model 1 under the hypothesis Hi that it is from model 1, 
while P(S2\Hi) is the probability to incorrectly identify a sample from model 2 under the same 
hypothesis. The probability to obtain the correct conclusion is thus 

n = (3) 

PiSilH^+PiSzlHi)- {> 
If the likelihood analysis of a CMB data set yields the conclusion that the data are drawn from 
model 1 and IT ~ 1, we can be confident we have made the correct conclusion. If IT ~ \ then no 
conclusions can be made concerning the result since we only have a 50-50 chance of being right. 



3. Models 

Any maximum likelihood analysis assumes the comparison of at least two competing models. 
We use three statistical tests (genus, extrema correlation, and bispectrum) of the COBE-DMR 
4-year sky maps to compare a Gaussian model against a semi-analytic texture model, a physical 
non-Gaussian theory. We generate 1500 simulated sky maps for each model, and use these maps 
to generate the mean values (S a ) and the covariance matrices 

M afS = 1 J>« - (S a ))(S^ - (S P )), (4) 

i 

including effects of instrument noise and Galaxy cut. We then generate an additional 1500 simulated 
sky maps of each of the models for a maximum-likelihood comparison, which we use to quantify 
the power of each statistical test. 
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The microwave anisotropy may conveniently be expressed as a sum of spherical harmonics, 
ATi = ^2i m aim Y[ m (9i, (pi). For the Gaussian model, we randomly distribute the spherical harmonic 
coefficients as a/ m = r exp(i5), where 5 is uniformly distributed between and 2ir and r a Gaussian 
random variable with zero mean and variance 

/i\-^L(n , 2 r(/ + (n-l)/2)r((9-n)/2) 

\r ) - 5 Wrms j r (/ + (5 _ n)/2) r ((3 + n)/2) • W 

We use the best-fit spectral index n = 1.2 from the COBE 4-year data (Bennett et. al. 1996; 
Gorski et. al. 1996; Hinshaw ei. al. 1996a; Wright et. al. 1996). Our results are insensitive to 
changes of 5n ± 0.2, roughly the uncertainty in the fitted index. 

Our semi-analytic model is based on probability distributions for the texture energy density. 
We use numerical simulations of the texture field equations to derive these distributions and in 
turn use these distributions to generate simulated CMB anisotropy sky maps. The main addition 
to previous semi-analytic models (Turok and Spergel 1990) is the inclusion of partial wound events 
(Borril et. al. 1994). In (Phillips and Kogut 1995), it was verified that these events, though less 
energetic than the unwinding events considered in (Turok and Spergel 1990), are more numerous. 
When they are included, the number of free parameters of the texture model increases from just 
the symmetry breaking scale to now include the minimum and maximum strengths of the partial 
wound events. As we discuss in the Appendix we take the maximum strength to be the same as 
the unwinding strength, since these are topological configurations that just barely become knots. 
The minimum strength is put at 1/4 the maximum. Our results are insensitive to the precise value 
chosen. 

We smooth each realization of a theory sky map with a Gaussian profile of 7° FWHM to account 
for the COBE-DMR beam profile. We add a pixel-wise noise to each map, T" lstr = ATi + n i, using 
the noise template appropriate for the 4-year COBE data as corrected for Galactic emission using 
either the "correlation" or "linear combination" models of (Hinshaw et. al. 1996b). All statistics 
are evaluated only for the high-latitude sky (3881 pixels at |6| > 20° with custom cutouts near 
Ophiuchus and Orion, (Bennett et. al. 1996; Banday et. al. 1996)) from which a fitted dipole and 
quadrupole have been removed. 



4. Genus 

The statistical properties of the CMB can be characterized by the excursion regions enclosed 
by isotemperature contours. The genus is the total curvature of the contours at fixed temperature, 
and may loosely be defined as the number of isolated high temperature regions (hot spots) minus 
the number of isolated low temperature regions (cold spots) (Gott et. al. 1990). The genus per 
unit area is a locally invariant quantity and is insensitive to incomplete sky coverage (e.g. removal 
of the Galactic plane). 

We use a nearest-neighbor search to count the number of isolated hot and cold spots separated 
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by temperature threshold v. For each map, we compute the genus at 31 values from v = —3a to 
v = +3(7, where a = \J (AT 2 ) is the standard deviation of the noisy map. These values form the 
vector S in the likelihood analysis. We adopt the difference in spot number as an estimate of the 
genus and use the Monte Carlo simulations to calibrate any difference induced by the Galactic cut 
between this definition and the total curvature. The two definitions are identical in the absence of 
a Galactic cut, and the nearest-neighbor algorithm allows significant computational savings. 

We use the genus to do a maximum likelihood fit for the model amplitude for the COBE 
data. We find that both models yield similar amplitudes, A^Gauss = 0-64 (Q r ms = 13.0 fiK) and 
^4fit,texture = 0.62 (Gi] 2 = 3.5 x 10 -7 ). The likelihood of the COBE data is greater for the Gaussian 
model. To derive the statistical power for the genus, we present the likelihood functions with 1500 
independent sampling maps, with amplitudes ranging from A = to A = 3 in 31 uniform steps. 
We fit each sample for the amplitude, then select the more likely model. We then filter these 
samples to have an amplitude close to the fitted DMR map: \Ast — ^4fit,Gauss| < 0.05 (The filter 
width corresponds to our sampling interval in A; our results for IT are insensitive to the exact 
value). We compare the number of input models identified as each model (Table 1). We obtain 
IlGauss = 990/(990 + 704) = 58%. Although the Gaussian model is the more likely model for the 
COBE data, we can not be overconfident about this result. When the noisier "combination" map is 
used, the statistical power drops to IlGauss = 51% and ITrexturc = 54%. Table 2 shows the statistical 
power for the various statistical tests. 



5. Extrema Correlation 

The extrema (peaks and valleys) of a field are defined as those points with VT = 0. For 
a pixelized map, this reduces to the collection of pixels hotter or colder than all of their nearest 
neighbors. Identifying pixels hotter than their neighbors produces a set of "peaks", while the locally 
coldest pixels produce the "valleys" . 

The number density of peaks or valleys is dominated by the noise properties of the map 
(Kogut et. al. 1995). The clustering of the extrema, as measured by the 2-point correlation function 
of the maxima and minima, provides additional information on the underlying CMB temperature 
field. We define the extrema correlation function as the 2-point function of all extrema points, 



Cext(0) 



Eij w i w j 

(Kogut et al. 1995), where the sum over pixel temperatures is restricted to the set of maxima and 
minima pixels separated by 9 evaluated on the cut sky. Since the correlation properties of the 
non-uniform noise in the DMR maps are different from the underlying CMB temperature field, 
we smooth the maps with an additional 7° FWHM Gaussian prior to collating the extrema as a 
compromise between suppressing noise and removing CMB power at small scales, as was pointed 
out in (Kogut et. al. 1995). 
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The clustering of extrema depends on the threshold u, as defined for the genus. We evaluate 
C ex t(#) at thresholds |z/| = 0, 1, and 2. Since, by definition, two peaks can not be adjacent, we 
ignore both the bin at zero angular separation and the first non-zero bin in all subsequent analysis. 
Simulations show that the results are dominated by the first few remaining bins; consequently, 
we speed processing by truncating the correlation function at separation = 60° for a total of 
20 angular bins of 3° per bin at each of the three thresholds. We concatenate the three binned 
extrema correction functions together to form our vector S for the likelihood analysis. A maximum- 
likelihood analysis of the COBE 4-year data shows the Gaussian model to be modestly more likely 
than the texture model. The statistical power of the extrema correlation is comparable to the genus 
(Table 2). 



6. Bispectrum 

A third way to characterize a random field is its hierarchy of n point functions. Since a Gaussian 
random field is completely characterized by its one and two point functions, the three point function 
or its harmonic conjugate, the bispectrum, may prove sensitive to non-Gaussian features in the data. 
We assume that the three point function only depends on the angular separation between points, 

C^(9 12 , 9 13 , 9 23 ) = (AT(m) AT(n 2 ) AT(n 3 )> , (6) 

where cos 6ij = • n,. The expansion in terms of Legendre polynomials (Gangui et. al. 1994) is 
given by 

h,h,h mi,m 2 ,m s 



j,k,lmj,m k ,mi 



with 



J Yf imi (n) Y hrn2 (n) Y hma (n) .dn (8) 

The bispectrum expansion coefficients are defined in terms of the spherical expansion coefficients 
aim via 

Our assumption that depends only on the angular separation Oij implies we are working 
with an isotropic cosmology. This leads us to consider the angular averaged bispectrum coefficients 

«.«.<.= E ('*. '! i)«KT w 



mi,m2,wi3 



mi 7TT-2 ITI3 
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where (• • •) are the Wigner 3-j symbols. We also introduce the normalized coefficients 

thhh - a h h h — - i~~t y 11 ) 



a hhh 



(C h C h C; 3 ) 2 

_1 ( h h h 

" '000 



(((2/i + 1)(2Z 2 + l)(2/ 3 + 



The diagonal (l\ = I2 = h = I) coefficients are the if coefficients of (Ferreira et. al. 1998) while 
the Jf coefficients of (Magueijo 2000) correspond to l\ = I — 2, I2 = I and I3 = I + 2. 

To compute the bispectrum coefficients for our likelihood analysis, we fit the spherical expan- 
sion coefficients a/ m on the cut sky up to a maximum value / max . We use these best fit a\ m to 
construct the bispectrum coefficients (10) and form the vector S = {lfu} l=2 4 ; • With this 
choice, our likelihood analysis uses the same statistic as (Ferreira et. al. 1998) and we find there 
is no statistical power: IIf3 jGauss = n_i-3 jtexture = 50%. Ferreira et. al (1998) used l max = 18; we 
find no change in statistical power as Z max is increased to 30. We find a similar lack of power using 
S = {|-B/h|}. We conclude that neither the bispectrum nor the normalized bispectrum in this form 
can differentiate between Gaussian and texture models, even though these models have statistically 
identifiable differences in the genus and extrema correlation. 

There are two possible explanations for this failure: i) the galaxy cut is aliasing too much 
power; or ii) we are not using enough of the available information present in the bispectrum. To 
address these possibilities, we vary both the size of the galaxy cut and the number of the off-diagonal 
terms of the bispectrum. For a statistical isotropic field, the two-point function has no coupling 
between different I multipoles: (a/ m a*/ m i) = C\ 5u> 5 mm i. We find this assumption can not be made 
for the three point function. Figure la shows the distribution of (l-Bzj z 2 z 3 |) values for an n = 1 
Gaussian model without noise or a galaxy cut (averaged over 1500 realizations). The diagonal 
values, (\Bui\), are dominant, but the off-diagonal terms also have non-zero values comparable to 
the diagonal. Moreover, the values group according to how much they are off the diagonal. There 
are three groups: the diagonal values, the terms off just one diagonal ((\Bhji^i\)), and finally the 
terms completely off the diagonal 1 ((|i?/ 1 1 2 1 3 1), ^ / I2 / I3). Figure lb shows the distribution 
for our texture model. The bispectrum coefficient values still cluster, but not as tightly as the 
Gaussian model. The Jf coefficients of (Magueijo 2000) are in this last group. Their choice of 
statistic misses the more significant set of terms off just one diagonal. The normalized bispectrum 
coefficients (^j 2 j 3 ) demonstrate similar clustering properties. 

The Galactic cut destroys the bispectrum's ability to differentiate models. Figures lc and Id 
show the distributions of (l-B^ j 2 <, 3 |) for Gaussian and texture models after a galaxy cut is included. 
The values no longer cluster, and are no longer identifiable by model. The spherical harmonics Y\ m 
are no longer orthogonal on the cut sky; the resulting aliasing of power between modes overwhelms 



1 We only consider the terms that are not forced to zero by the triangularity of the 3-j symbols. 
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the signal. 

To quantify this failure, we compute the statistical power as we impose the simple latitudinal 
cut of removing pixels with 9 < b, for b varying from 0° to 14°. We set a signal amplitude 
that corresponds to the genus and extrema correlation function (A = 0.64). Since the bispectrum 
computed using only the diagonal values never shows statistical power, we include both the diagonal 
and once-off-the-diagonal terms (the two rightmost groups in Fig 1). Figure 2 shows that the 
probability to correctly identify the input model falls monotonically as the Galactic cut is increased. 
For comparison purposes, we also plot the statistical power for the genus, which shows only a slight 
trend consistent with larger sample variance as the solid angle available for analysis is reduced. 
With no galactic cut, the bispectrum has comparable power as the genus 2 . The statistical power 
of the bispectrum drops steadily until cut b < 12°, where there is no power. The normalized 
bispectrum 1^ v is never able to differentiate between the models, even when there is no galaxy cut 
present. 

The other important consideration is how much information is needed. For / max = 20, there 
are 770 nonzero, unique bispectrum values, of which 10 are the non-zero diagonal values. Using all 
values instead of just the diagonal values is computationally expensive. We use instead the diagonal 
terms and terms with one off-diagonal I value (second group from the right in Figure 1) to form 
the vector S = The restriction \l — l'\ < Al selects how many of the off-diagonal terms 

we use. Al = corresponds to the analysis of (Ferreira et. al. 1998) and Al = l max includes all the 
once-off-the-diagonal terms. Figure 3 shows the statistical power of the bispectrum as we increase 
Al, where we consider both noisy and noiseless data sets evaluated over the full sky. The power 
is greater without noise; the statistical power of the Gaussian model is fairly constant after the 
inclusion of just a few extra off-diagonal terms. For the texture mode, the inclusion of additional 
off-diagonal terms steadily improves the probability for correct identification until Al ~ Z m ax/2. 
With noise, the statistical power for either model levels out at about the same value as the genus 
and extrema correlation function. This suggests that but for the need to include a galaxy cut, a 
likelihood analysis based on \Bm>\ would be as powerful as the other tests. The galaxy cut aliases 
power in the conjugate harmonic space, but has negligible effect in the coordinate space used by 
the genus and extrema correlation function. 

The addition of off-diagonal terms is important for the non-Gaussian texture model. For 
ideal (noiseless) data, the probability to correctly identify the texture model goes from 59% for 
only diagonal terms to 74% when all \Bmi\ terms are considered. In contrast, the probability to 
correctly identify the Gaussian model is nearly constant. Since the three point function contains 
no new information for a Gaussian theory, we do not expect much improvement in the Gaussian 
model as the number of terms is increased. The texture model was chosen exactly because it is an 
example of a non-Gaussian theory, so the three point function or bispectrum should be sensitive to 
the gradual addition of the information present in the off-diagonal terms. 



2 Since we fix the input amplitude instead of fitting, the genus power here is slightly larger than shown in Table 2. 
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We have duplicated all the analysis discussed above for the normalized bispectrum coefficients 
introduced by (Ferreira et. al. 1998). Figs 2 and 3 show the likelihood analysis in terms of 
the normalized bispectrum never has any statistical power. For the ideal situation of noiseless 
data without a galaxy cut, the probability to correctly identify either model is at best 52%. The 
statistic Jf of (Magueijo 2000) corresponds to our lf_ 2 1 1+2 and is thus not directly covered in 
our analysis. Figures la and lb show jf to be of lesser significance than our most general set of 
ifur's and we can not expect any statistical power for this choice. Although we have tested only a 
single non-Gaussian model, the normalized bispectrum does not appear to be a promising test of 
the statistical distribution of CMB anisotropy. 



7. Discussion. 

Recent work has generated new interest in statistical tests of the CMB. Comparison of the 
bispectrum of the COBE-DMR 4-year sky maps to simulations of Gaussian maps show the COBE 
data to lie far from the mode of the Gaussian models (Ferreira et. al. 1998). Since this work tested 
only Gaussian models, it is not clear whether the discrepant values result from non-Gaussian signals 
in the data (including instrumental artifacts) or are simply an outlier drawn from a Gaussian parent 
population. We compute the statistical power of several popular statistics (the genus, extrema 
correlation, and bispectrum) to determine their ability to distinguish between a Gaussian model 
and a texture model chosen as a physically motivated non-Gaussian alternative. 

Both the genus and extrema correlation successfully discriminate between the models with 60% 
confidence, limited by the similarity of the texture model to a Gaussian on large angular scales (the 
central limit theorem). The COBE 4-year data prefer the Gaussian model, although the modest 
statistical power prevents a strong rejection of this specific non-Gaussian alternative. In contrast, 
we find that the bispectrum has no power to discriminate between models, and trace this to the 
aliasing of power caused by the Galactic cut imposed on the data. The normalized bispectrum 
never has any power, even run on noiseless simulations with no Galactic cut. Since the results 
presented in (Ferreira et. al. 1998) and (Magueijo 2000) are based on an analysis in terms of the 
normalized bispectrum, the results, though correct, have no statistical power. Their results provide 
a necessary, but not sufficient, condition to rule out an underlying Gaussian character to the CMB 
anisotropy. 

The problem of aliased power can in principle be solved, either by windowing the multipole 
coverage or by construction of an explicitly orthogonal basis. With this in mind, we further explore 
the statistical power of the bispectrum using full sky coverage to avoid the problem of aliased power. 
We find significant coupling between multipoles in the bispectrum coefficients for the non-Gaussian 
(texture) model. Failure to include these off-diagonal terms reduces the statistical power of the 
bispectrum on the full sky from 67% (comparable to the genus or extrema correlation) to 55% (not 
much better than random guessing). The off-diagonal elements are more important for the non- 
Gaussian model, pointing out the importance of testing various statistics against some alternative 



to the standard Gaussian model. 
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Appendix 



A. 



Texture Model 



In order to test the power of different statistics, we need an alternative hypothesis to the stan- 
dard Gaussian model. Defect models (Kibble 1985a; Kibble 1985b; Vilenkin and Shellard 1995) 
are the only known cosmological alternatives expected to provide a significant non-Gaussian signal. 
Characterizing this non-Gaussian signal requires many realizations of the model. Unfortunately, 
the non-linear nature of defect models make it impossible to carry out the analysis in closed form, 
while also making exact numerical work computationally expensive. We use exact simulations of 
a texture model to derive statistical descriptions of the energy density, then use these statistical 
descriptions to rapidly generate additional realizations of the model. 

Textures (Turok 1989) are a class of topological defects in which the vacuum manifold of the 
order field in broken symmetric phase is of the same dimension as the spatial geometry. This allows 
the order field to remain on the vacuum manifold, regardless of initial conditions. Perturbations in 
the energy density are driven solely by the order field's kinetic energy. As the correlation length 
of the order field grows, causal regions emerge that posses non-trivial topological charge (knots). 
These configurations collapse until the energy density is great enough for the order field to tunnel 
through the vacuum manifold and unwind. As CMB photons travel through these evolving regions 
of increased energy density, they become blue/red shifted and the anisotropy arises. This basic 
picture of knots in the texture order field led to the first approximation of Turok and Spergel 
(1990) which randomly placed idealized knots across the sky. The resulting CMB anisotropy, 
however, did not agree with results derived from full numerical simulations (Pen et. al. 1994). 

Borril et. al. (1994) suggested the existence of a less energetic but more numerous configura- 
tion: partially wound events (PWEs), in which order field wraps around the vacuum manifold but 
not enough to be knotted. In (Phillips and Kogut 1995), we directly verified, for an expanding flat 
universe, that PWEs are far more numerous than knots and that their contribution dominates the 
CMB anisotropy. In this Appendix, we provide a statistical description of the PWEs. We then ex- 
tend the model of (Turok and Spergel 1990) by randomly placing texture events on the microwave 
sky, drawing each event from the correct distribution in energy. 



Since we are interested in times long after the phase transition for the order field, we treat the 
field as classical. Except for the special case where the field unwinds, we restrict the field to lie on 
the vacuum manifold. We will focus solely on textures, so the (real-valued) global order becomes 
4? = (4>i, $2, $3, $4 ) and the vacuum manifold is the three sphere S* 3 : |4?| 2 = 4>q. After re-scaling 



A.l. Field Equations and Numerical Implementation 
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the field <fr = &/4>o, the field dynamics are determined by the action 

"1 



m-J 



d Xy/—g 



(a M *)-(^*)-A(|*| 2 -i) 



(Al) 



The constraint that the field stay on the vacuum manifold is imposed via the Lagrange multiplier 
A. This is the nonlinear a model (NLSM) (see e.g. (Pen et. al. 1994)). We carry out our work for 
a spatial flat Robertson- Walker homogeneous universe. For this choice the field equations become 



* + 2-* - V 2 * 

a 



(A2) 



where <& = d<&/dr and V3> = d^/dxi is the spatial gradient. We refer the reader to (Phillips and Kogut 1995) 
for details of the derivation of discrete versions of these equations and their numerical implemen- 
tation. 

We characterize the field dynamics using two quantities. The first is the energy density, the 
tt component of the field stress tensor: 



Pi,n 



1 



i.ii 



+ |V*i,J 



(A3) 



with the obvious discrete versions of the derivatives, where index i labels the spatial grid cells and 
t is the temporal label. The second quantity is the cell alignment, 



(A4) 



where j runs over cell i's six nearest neighbors. With this definition, the alignment of a grid cell 
is the cosine of the angle between the field configuration at i and the average of the field over i's 
nearest neighbors, aj measures how close the vacuum at % is to its neighbors, or how aligned the 
cell is with the local choice of vacuum. This gives a rough measure of the local topological charge 
(see (Borril 1995) for further discussion). According to this use of the NLSM model, knots are 
numerically identified by on < (Pen et. al. 1994). They are unwound by flipping the field to its 
anti-podal: 3>j — > — 

We base our analysis on a catalogue of events identified in 4000 simulations of the order field 
dynamics. The expansion scale factor a(t) corresponds to a matter-dominated universe. We run 
the simulations on a 64 3 grid array with periodic boundary conditions, continuing until the horizon 
crossing time (At = 32). The initial order field configuration <£(x, To) is uniformly distributed over 
the vacuum manifold. We begin the simulation by integrating this initial random configuration, 
but with the anti-podal flipping of knot configurations disabled. This acts to "smooth" the initial 
order field while maintaining the constraint |3>| 2 = 1. We then re-enable the anti-podal unwinding 
and begin the event identification procedure outlined below. 
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A.2. Event Identification 

As the casual horizon grows and regions of different vacua come into contact, the order field 
$(x) changes until it is in the same vacuum throughout. As this takes place the energy density- 
due to the field rises. The dynamics of <1> force the region containing the changing $ to shrink 
in size (Derrick 1964). When the total topological charge contained in this region is non-trivial, 
i.e., a knot, the energy density within the shrinking region will increase until the field can tunnel 
through the vacuum manifold and unwind the topological charge. The other possibility is that the 
order field within the region, though needing to resolve which vacuum the region will occupy, has 
no net topological charge. These are the partially wound events. They too are shrinking regions of 
growing energy, but now order field sorts itself out without having to tunnel through the vacuum 
manifold. In either case, after the region shrinks and resolves which vacuum to occupy, the region 
then expands as the stress energy diminishes. 

We use this behavior to identify events in the numerical simulations. At each time step, we 
record the position, alignment, and energy density for each local maximum in the energy density, 
then track the spatial movement of each local maximum from one time step to the next. These 
local maxima are the grid centers of the growing/shrinking regions discussed above. We identify 
300597 events in the 4000 simulations run. Figure 4 shows the temporal evolution of typical events. 
The alignment decreases and the energy increases as the texture becomes increasingly wound. 
Shortly after the alignment reaches a minimum, the energy density begins to decrease and the 
event dissipates. 

Fig. 5 shows the distribution of these events sorted by the niininiuni alignment c^min 

of each 

event. Three groupings are evident: knots, PWEs and noise. The knots are all events with minimum 
alignment less than zero. We identify all events with a m i n > 0.55 as noise, since they have negligible 
contribution to the energy of the field. We identify all remaining events, < a m \ n < 0.55, as 
partially-wound events. With this identification, 4% of the events are knots, 56% are PWEs, and 
40% are noise. 

A. 3. Partially Wound Event Distributions. 

A semi-analytic model for textures requires a statistical description of the temporal and energy 
distribution for the PWEs. Instead of working directly with the energy distribution, we consider 
the alignment of the PWEs. The alignment is more closely related to the topological properties of 
the theory and as such is less sensitive to noise in the simulated order field. 

We label each event by the minimum of its alignment, a m j n , and the conformal time r m i n of 
this minimum alignment. Binning the PWEs according to a m i n and r m i n shows the distribution to 
be separable, d/VpwE^, t) oc dNi (r) dAT 2 (a) . Not surprisingly, the temporal distribution for the 
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PWEs has the same behavior as that for the knots, 



diVi(r) = KpwE dr; k PW e = 7- (A5) 

5 



We determine kpwe by comparing the total number of PWEs (iVpwE = 168307) to the number 
knots (iVknot = 11295). There are 15 partially wound events for each unwinding event. 

Figure 6 shows the distribution in alignment, diV^cc). The shape is well described by a x 2 
distribution with 4 degrees of freedom, 

dN 2 (a) = A[xl(x)+p ] (A6) 

where xl( x ) ls the x 2 distribution for v = 4 degrees of freedom, 

a — a 

x = , A7 

m 

a = 0.55 is the alignment noise cutoff, and A is a normalization constant. The uncertainty in the 
number of events Nj in each bin is given by counting statistics, a 2 N . = Nj. We thus fit the binned 
data to determine values 



m = 0.165 ±0.003 and Po = 0.026 ± 0.005 (A8) 

A four-parameter fit with v and a free yields fitted values u = 4.01 ± 0.02 and a = 0.54 ± 0.01, 
justifying the choice of functional form of the probability distribution. 

We also test the spatial properties of the partially-wound events, to determine whether one 
PWE will suppress or enhance the probability to find another PWE in the local neighborhood. The 
two-point correlation function for all PWEs is indistinguishable from one derived for the same num- 
ber of events uniformly distributed in space. A principal moment analysis of the spatial distribution 
of the energy density around the gird center shows the PWEs to be spherically symmetric. 



A. 4. Sky Map Generation 

Our semi-analytic model uses this statistical description of the order field in texture cosmolo- 
gies to generate maps of simulated CMB anisotropy without requiring computationally expensive 
numerical simulation of the evolving order field. It extends the work of Turok and Spergel (1990) 
to include the more numerous partially- wound events. As in that work, we generate a random set 
of of spherically symmetric events with uniform spatial distribution and r -4 temporal distribution. 
Now, however, the energy of each event is drawn either from a mono-energetic knot population (6% 
probability) or from the distribution of PWEs (94% probability). 

Deriving a map of the CMB anisotropy from a realization of random texture events requires 
fixing the relation between the minimum alignment a used to generate the PWE distribution (Eq. 
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A6) and the amplitude AT/T of the CMB anisotropy. Analysis of the simulations indicates this 
relation to be linear. Due to the discrete anti-nodal flipping algorithm, the energy scale of the knots 
in the simulations is arbitrary. We thus require two normalization parameters to fully specify CMB 
anisotropy given a distribution of texture events: the minimum and maximum energy of the PWEs 
relative to the knots. We take the maximum amplitude for PWEs to be equal to the amplitude 
for the knots: these represent (rare) events that just barely missed unwinding (a ~ 0). The noise 
cutoff a m j n = 0.55 provides a lower limit to the energy of PWEs used in the model; we thus take 
the minimum amplitude to be half way between the knot amplitude and zero. 

A single realization of a CMB map using the semi-analytic model thus requires the following 
steps: 

1) we compute the number of events in the observed universe by taking the product of the 
observed volume and the integrated temporal distribution (A5); 

2) we distribute these events uniformly in space and with a r -4 temporal distribution; 

3) letting 6% of these be knots, we assign them the maximum amplitude and let the other 94% 
be PWEs. For the PWEs, we select random alignments from the distribution (A6) and linearly 
map to amplitudes as outlined above; 

4) we use the result from (Turok and Spergel 1990) for the AT/T contribution to determine 
which pixels are influenced by each event. 

We are thus left with one undetermined parameter, the amplitude of the knotted configurations. 
A single full-sky realization of the semi-analytic model with 6144 pixels requires only 20 seconds 
of CPU time on a Sun Ultra workstation at 360 MHz clock speed, achieving the goal of rapidly 
generating a large number of sky maps. 

We thank Ue-Li Pen for valuable discussions. This work was funded in part through NASA 
RTOP 399-20-61-01. 
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Fig. 1. — Distribution of average bispectrum values (|-Bj 1 j 2 j 3 |) for 1500 maps, (a) Noiseless n = 1 
Gaussian model, no galaxy cut. (b) Noiseless -E m i n = 0.25 E max texture model, no galaxy cut. (c) 
Same as (a), but now with a galaxy cut. (d) Same as (b), but now with a galaxy cut. (For ease of 
viewing, each of the bispectrum was rescaled to a common range. All analysis was done without 
this rescaling.) 

Fig. 2. — Dependence of the statistical power of the bispectrum l-B^^^I and power-spectrum 
normalized version ^ t on the size b of the galactic cut. The statistical power is measured by 
the probability to correctly identify the input model. 

Fig. 3. — Statistical power of the bispectrum \Bm'\ and normalized bispectrum I?,,, as a func- 
tion of the number of off-diagonal terms. The number of off-diagonal terms is controlled by the 
restriction \l — l'\ < Al. Al = uses ony the diagonal terms and is appropriate only if there is no 
coupling between multipole terms in the bispectrum. The texture model (triangles) has significant 
contribution from off-diagonal terms. The normalized bispectrum has no power regardless of the 
number of off-diagonal terms. 

Fig. 4. — Temporal evolution of typical texture events. The top panel shows the alignment at the 
center of each event at each time step. The bottom panel shows the energy density at the center, 
defined as the local maximum of the energy density. 

Fig. 5. — Distribution of texture events sorted by the minimum alignment reached over the history of 
each event. Three peaks are evident. Events with a m i n < correspond to knots. The energy density 
at the time of minimum alignment (bottom panel) shows that the region a m j n > 0.55 contributes 
negligible. We identify these events as noise. The middle region corresponds to partially-wound 
events (PWEs). 

Fig. 6. — The distribution function dN2(ce)Aa for partially-wound events. The smooth curve is 
the fitted xt distribution with 4 degrees of freedom. 
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Table 1. Distribution of Model Identifications 
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Input Model 


Gaussian 


Texture 


Gaussian 


990 


451 


Texture 
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845 
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Table 2. Statistical Power of Analyzed Tests 
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Noise Map 


Amplitude 


Galaxy Cut Gaussian 
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Corr 


As* = 0.64 


Template 
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Comb 


A fit = 0.66 


Template 
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A fit = 0.61 


Template 


61 
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A fit = 0.72 
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54 
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Bispectrum with varying Galaxy Cut 
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A = 0.64 
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None 
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None 
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None 
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51 


Bispectrum with and without Off-Diagonal Terms, with Noise 
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64 
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50 


50 
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51 


51 


Bispectrum with and without Off-Diagonal Terms, Noiseless 


Bm 
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68 


59 
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71 


74 


4i 
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50 
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52 


52 



Distribution dN( a )Aa 
o ^ ro co ^ 




